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SUMMARY 


A  comparison  is  made  between  experiment  and  theory  for  the  calcu¬ 
lation  of  shielding  by  a  structure  against  fallout  radiation. 

Experimental  results  for  exposure  penetrating  a  roof  slab,  and 
for  reduction  factor,  are  found  to  agree  with  moments  theory  calcu¬ 
lations,  often  to  better  than  10%  when  the  geometry  factor  L  (u),X) 

a 

was  used. 

A  comparison  between  experiment  and  theory  may  be  inaccurate  due 
to  anisotropy  of  the  experimental  source,  to  lack  of  source  reflection 
in  roof  experiments,  and  due  to  error  in  estimating  the  thickness  of 
the  roof.  The  magnitudes  of  the  errors  due  to  these  effects  are  inves¬ 
tigated,  and  found  to  be  small,  but  not  necessarily  negligible. 

Detailed  results  of  penetration  in  iron  and  concrete,  due  to  plane 
isotropic  sources  of  various  energies  are  given,  and  a  review  of  ground 
contribution  to  fallout  radiation  penetration  through  simple  structures 
is  appended. 
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1  INTRODUCTION 


Standard  methods  of  calculation  of  the  protection  afforded  by  a 
structure  against  fallout  radiation  depend  upon  the  separation  of  the 
solution  into  a  barrier  and  a  geometry  factor.  These  factors  have  been 

calculated  in  large  part  by  moments  theory,  and  have  been  presented  as 

.  (1,2) 

graphs 


The  primary  aim  of  the  research  project  reported  herein  was  to 
determine  how  valid  the  method  of  calculation  is  when  applied  to  the  case 
of  fallout  radiation  on  the  roof  of  a  structure.  In  order  to  do  this, 
experimental  work  on  the  subject  was  reviewed  and  compared  with  theory. 

A  number  of  reported  experiments  were  examined  and  it  was  found  that  much 
of  this  work  had  been  on  structures  such  as  houses  and  barracks  which  had 
very  complicated  geometrical  configurations.  Work  reported  on  such  struc¬ 
tures  gave  values  of  radiation  exposure  valid  only  for  the  structure 
tested  and  accurate  generalization  of  the  results  to  other  structures 
appeared  to  be  extremely  difficult.  Since  the  aim  of  this  project  was 
to  best  the  method  of  calculation  of  protection  against  fallout  radiation, 
only  those  experiments  which  had  been  conducted  on  structures  of  reason¬ 
ably  simple  geometry  were  used  to  compare  with  theory.  The  only  such 
experiment  reported  at  the  time  of  writing  was  that  at  the  Nuclear  Defense 

(3) 

Laboratory  .  Preliminary  results  from  work  at  the  Protective  Structures 

(4) 

Development  Center  were  also  available. 


To  provide  a  comparison  between  theory  and  experiment,  it  was  neces¬ 
sary  to  compute  the  exposures  due  to  a  plane  isotropic  source  penetrating 
steel  slabs,  since  steel  slabs  were  used  in  the  NDL  experiment  and  cal¬ 
culations  for  steel  using  moments  theory  were  not  available.  At  the  same 
time,  calculations  were  done  for  concrete  and  for  air,  for  a  plane  isotropic 
source,  at  several  energies. 

Three  factors  which  may  cause  error  in  comparison  of  theory  and  experi¬ 
ment  were  considered.  These  are  the  effect  of  source  anisotropy  in  free 
field  measurements,  the  effect  of  lack  of  source  reflection  in  roof  experi¬ 
ments,  and  the  effect  of  error  in  estimating  the  thickness  of  the  slab 
used  in  the  experiment.  The  order  of  magnitute  of  errors  due  to  the  above 
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three  factors  were  estimated  and  they  were  found  generally  to  be  small 
but  not  necessarily  negligible,  as  reported  in  Section  3. 

(3  4) 

The  experimental  measurements  available  ’  were  not  corrected  for 
these  effects  in  this  report.  The  experimental  results  were  compared 
with  the  theoretical  exposures  and  reduction  factors  calculated  here  (see 
Appendix  A);  agreement  was  often  better  than  10%. 

Tn  the  appendices,  two  subjects  of  related  interest  which  arose 
during  this  investigation  are  reported.  These  are  the  detailed  results 
of  exposure  distributions  and  reduction  factor  calculations  in  iron  and 
concrete  due  to  plane  isotropic  sources  of  various  energies,  calculated 
by  moments  theory;  and  a  review  of  ground  contribution  to  fallout  radiation 
penetration  through  simple  structures. 

2  THE  LOGIC  OF  THE  METHOD  OF  ANALYZING  STRUCTURES  FOR  FALLOUT  PROTECTION 

Before  proceeding  with  a  comparison  of  theory  and  experiment,  it  is 
helpful  to  know  where  the  theory  fits  into  engineering  technology  and 
how  experiment  relates  to  the  theory.  The  relation  of  theory  and  experi¬ 
ment  are  discussed  in  this  Section,  together  with  some  general  considerations. 

In  engineering  design,  the  properties  of  a  structure  or  other  item 
are  specified,  and  the  components  are  then  selected  and  proportioned  so 
as  to  give  the  required  results;  in  analysis,  a  structure  is  completely 
specified  and  then  some  property,  such  as  its  behavior  under  an  unusual 
condition,  is  calculated. 

If  a  calculation  method  is  available  which  permits  direct  computation 
of  the  size  and  specifications  of  the  structure  to  satisfy  the  design 
criteria,  then  the  structure  may  be  directly  designed.  Many  structures 
cannot  be  directly  designed.  The  procedure  then  is  to  select  a  design, 
analyze  it,  amend  the  design  to  satisfy  some  given  criteria,  analyze  again, 
and  so  on.  We  then  have  a  feed-back  interative  process  where  analysis 
is  one  part  of  the  design  procedure. 

The  hypothetically  most  accurate  and  direct  approach  to  the  problem 
of  analyzing  a  structure  for  its  ability  to  shield  people  from  fallout 
radiation  would  be  to  subject  it  to  actual  fallout  and  measure  the 
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exposure  received.  This  Is  not  practical,  nor  is  it  practical  to  subject 
every  structure  to  simulated  fallout  radiation. 

For  this  reason,  structures  are  divided  into  classes,  categorized 
by  some  general  properties  such  as  size,  construction  material,  etc. 

This  implies  that  we  know  which  properties  are  significant  for  this  prob- 
4  lem,  and  which  are  not,  so  that  the  process  of  generalization  has  small 

effect  on  the  accuracy  of  the  analysis. 

It  is  not  practical  to  do  experiments  which  cover  all  relevant  regions 
of  all  relevant  variables  even  within  one  class  of  structure.  Recourse 
is  then  made  to  theory,  which  is  a  generalization  of  the  practical  problem 
in  a  way  that  can  be  applied  with  reasonable  effort  and  accuracy  to  any 
structure  in  the  class  to  which  it  pertains. 

Before  the  theory  of  radiation  protection  can  justifiably  be  used, 
it  must  be  established  for  which  situations  it  is  applicable  and  how 
accurate  it  is.  This  is  done  by  comparing  theory  with  experiment,  and 
,  the  experiment  must  be  as  careful  and  as  accurate  as  possible  so  that  it 
gives  some  evidence  as  to  the  validity  of  the  theory.  Once  that  is 
established,  the  theory  may  be  used  with  an  established  degree  of  confi¬ 
dence,  for  any  particular  problem. 

Since  the  biological  effects  of  radiation  dose  are  relatively  insen¬ 
sitive  to  small  changes  of  dose,  great  accuracy  in  the  knowledge  of  the 
shielding  capability  of  any  one  structure  is  superfluous.  If,  therefore, 
an  experiment  were  performed  on  a  structure  simply  to  obtain  the  protec¬ 
tion  factor  of  that  structure,  it  need  not  be  extremely  accurate.  But 
if  it  is  performed  to  check  a  theory,  it  must  be  well  understood  and 
carefully  carried  out,  so  as  to  give  limits  of  confidence  to  the  theory. 

By  this  means  a  small  number  of  careful  experiments,  linked  to  a  theory, 
can  replace  much  random  empiricism. 

In  order  to  draw  conclusions  about  the  method  of  calculation  for 
analyzing  structures  against  fallout  radiation,  it  is  convenient  to  con¬ 
sider  both  the  experimental  and  practical  situations,  and  the  theoretical 
model,  in  five  parts  (see  Table  1): 
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(lit  The  spatial  distribution  of  the  source; 

(2)  The  uuciear  properties  ol  the  source; 

.3)  The  spatial  distribution  of  the  structure; 

(4)  The  nuclear  properties  of  the  structure; 

(5)  The  definition  of  the  dependent  variable,  exposure. 

These  five  aspects  are  discussed  below. 

(1)  It  is  assumed  that  real  fallout  will  be  uniformly  distributed 

in  plan.  The  validity  of  this  assumption  is  not  discussed  in  this  report. 
The  source  is  treated  as  a  uniform  plane  in  theory,  and  this  can  be 
accurately  simuiaitd  by  experiment. 

(2)  The  interactions  of  gamma  radiation  with  matter  are  well  under¬ 
stood  over  the  relevant  energy  range.  Conclusions  drawn  about  the  theory 
on  the  basis  of  experiments  with  sources  other  than  fallout  may  therefore 
be  taken  as  valid  for  application  to  real  situations.  There  is  some  dif¬ 
ficulty  in  knowing  just  what  to  consider  to  be  the  source  properties  in 
an  experiment,  on  account  of  the  fact  that  experimental  sources  are  often 
large  relative  to  the  mean-f ree-path  length  of  radiation  within  them,  and 
because  in  calculations  using  moment  theory  to  find  the  penetration  through 
a  roof  slab  the  exposure  originating  from  fallout  on  the  ground  (skyshine) 
is  not  separated  from  the  exposure  originating  from  fallout  on  the  roof. 
This  is  discussed  in  Section  3. 

(3)  The  theory  assumes  the  roof  of  the  structure  to  be  an  infinite 
homogeneous  slab  whereas  real  structures  consist  of  columns,  walls,  beams, 
and  finite  slabs.  'his  is  the  major  point  of  different e  between  theory 
and  reality,  and  experiments  should  be,  and  are,  designed  to  check  the 
validity  of  the  theory  in  different  geometries  of  structures. 

The  material  geometries  of  the  experiment  correspond  to  those  of 
real  structures.  In  order  to  use  experimental  evidence  to  draw  conclu¬ 
sions  about  the  relevance  of  theory  to  reality,  differences  in  source 
spectra  between  experiment  and  theory  must  be  accounted  for. 

(4)  Nuclear  properties  of  structural  materials  with  respect  to  ganuna 
radiation  are  known.  The  effect  of  the  small  Inaccuracies  in  available 
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nuclear  data  is  much  less  than  the  effect  due  to  the  relatively  large 
dimensional  tolerances  allowed  in  the  sizes  of  structural  members. 

(5)  The  detectors  used  in  experiment  measure  exposure.  This  is  the 
dependent  variable  in  the  theory.  It  is  worth  pointing  out  that  the 
absolute  calibration  error  of  even  the  best  detector  is  of  the  order  of 
2%  or  more,  and  some  detectors  used  have  inaccuracies  an  order  of  magnitude 
greater. 

Table  1  shows  in  summary  form  the  areas  of  similarity  and  dissimi¬ 
larity  as  discussed  above. 


TABLE  1 

Summary  of  relation  of  theory  and  experiment  to  assumed  real  situation. 
0  represents  similarity  to  assumed  real  situation 
X  represents  difference  compared  with  assumed  real  situation 


Property 

Theory 

Experiment 

(1)  Source,  Spatial 

0 

\ 

X 

(2)  Source,  Nuclear 

0 

X 

(3)  Structure,  Spatial 

X 

0 

(4)  Structure,  Nuclear 

0 

0 

(5)  Detector  response 

0 

0 

3  ANALYSIS  OF  SOME  POSSIBLE  CAUSES  OF  EXPERIMENTAL  ERROR 
3.1  Sources  and  Source  Calibration 

The  reduction  factor  of  a  structure  is  expressed  as  the  ratio  of  the 
exposure  observed  at  a  given  point  in  a  structure  to  the  exposure  which 
would  be  observed  three  feet  above  a  smooth,  infinite,  unobstructed, 
uniformly  contaminated  plane.  Experimentally  these  exposures  are  derived 
from  separate  measurements.  Each  exposure  is  normalized  to  an  arbitrary 
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source  strength,  usually  1  curie/ft  .  and  the  ratio  then  taken.  The 
protection  factor  is  defined  as  the  reciprocal  of  the  reduction  factor. 

If  the  same  source  were  used  for  the  structure  and  the  f  .-ee  field 
expe'iment,  and  if  it  were  isotropic  there  would  be  no  need  to  normalize 
the  dose  to  a  source  strength,  because  the  exposure  observed  in  the  struc¬ 
ture  divided  by  the  free  field  exposure  would  be  the  protection  factor. 

I f  an  anisotropic  source  were  used  in  either  exp2riment,  normalization 
of  the  observed  exposure  to  a  calibrated  source  strength  would  require 
that  the  source  anisotropy  be  allowed  for. 

In  practice,  the  sizes  of  the  sources  and  capsules  used  are  not 
small  compared  with  the  mean- f ree-path  length  of  radiation  in  them,  nor 
are  they  spherically  symmetrical.  This  leads  to  two  effects: 

(a)  Degradation  of  the  source  spectrum: 

(b>  Anisotropy  of  the  emitted  exposure  with  respect  to  the  source 
orientation. 

When  the  source  calibrated  in  air,  the  low  energy  components  of 
the  spectrum  are  detected.  When  the  source  is  placed  on  a  barrier,  the 
low  energy  components  are  more  rapidly  attenuated  than  the  higher  energy 
components,  If  the  results  are  compared  with  a  theory  derived  for  the 
source  energy,  the  theory  will  overestimate  the  transmitted  exposure. 

( 5 

Berger  and  Doggett  compared  thetr  calculations  with  an  experi- 

(61 

ment  of  Kirn,  Kennedy  and  Wyckoff  and  corcluded  that  neglect  of  the 
energy  degradation  from  a  1  cm  thick  Co-60  source  gave  rise  to  an 
ovemalculation  in  exposuie  of  18  .  for  barriers  approximately 

120  psf  thick. 

The  error  introduced  in  an  experimen*  by  neglecting  the  degraded 

source  spectrum  can  be  estimated  ll  the  source  spectrum  is  known,  by 

using  the  penetration  calculations  for  various  incident  energies  given 

in  this  report  in  Appendix  A.  Unfortunately  the  source  spectra  were 

not  determined  lor  sources  used  in  the  experiments.  This  effect  may 

tn  fact  turn  out  to  be  small  in  some  experiments  but  no  firm  evidence 

« 

is  available  to  allow  this  assumption. 

* 

A  program  to  calculate  by  Monte  Carlo  methods  the  spectr  as  of  gamma 
radiation  emitted  from  a  source  is  at  the  present  time  being  written, 
and  when  available  l*  should  help  to  solve  this  problem. 
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The  experimental  value  of  exposure  in  air  three  feet  above  the  source 
plane  is  obtained  by  measuring  the  exposure,  either  when  the  source  is 
pumped  through  a  tube  on  the  ground,  or  by  placing  a  detector  in  many 
locations  when  a  source  is  placed  in  one  location  and  the  results  summed 
to  Simula*  a  plcne  field.  An  estimate  of  the  effect  oi  source  anisotropy 
in  the  latter  procedure  is  made  in  Section  3.2. 

3.2  Effect  of  Source  Anisotropy  in  Free  Field  Measurements 

(7) 

Th'  free  field  measurements  made  at  the  Nuclear  Defense  Laboratory 
had  been  carried  out  by  placing  detectors  in  a  number  of  positions  relative 
to  a  single  source  on  a  field  and  summing  the  observed  exposures  to  s-nru- 
latt:  an  infinite  plane  source.  Since  the  source  was  probably  not  iso¬ 
tropic,  the  observed  exposure,  when  normalized  to  the  calibrated  source 
strength,  may  have  been  in  error.  An  estimate  of  the  error  due  to  this 
effect  was  made  as  follows. 

Since  the  anisotropy  of  the  source  radiation  was  not  known  for  the 
-  ources  used  at  the  Nuclear  Defense  Laboratory,  an  estimate  of  the  order 
of  magnitude  of  the  effect  was  made  for  a  hypothetical  case,  using  the 
reported  exposure  distribution  around  a  Co-60  source,  given  by  Burson 
and  Summers'  and  shown  in  Fig.  1. 

The  reduction  factor  due  to  area  n  of  Co-60  in  a  field  is  given  by 


~  -  P (d)  A/4rrd2  3.2.1 

o 


where  A  and  d  are  shown  in  Fig.  2  and  P(d)  for  Co-60  is  given  in 
Fig.  B-19  of  the  monograph  by  Spencer  Referring  to  Fig.  2,  this  can 

be  written  for  a  plane  isotropic  source 


D 

D 


P(d)  2nr  dr/4nd“ 


3.2.2 


The  exposure  emitted  by  the  Co-60  source  shown  in  Fig.  1  can  be  written, 
with  some  approximation  as 
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J  i 1  -  0.0033  G  )  . 

[  °-5 
Including  this  anisotropy  effect  in  eq. 


-60° 

<  9^60° 

3. 2.3a 

60° 

<  9^90° 

3.2.3b 

3.2.2  gives  an  approximate  relation 


D 

D 


F(6>  P(d)  r  dr/2d2 


3.2.4 


where  D  is  now  the  exposure  at  the  detector  due  to  a  uniform  infinite 
contamination  corresponding  to  a  calibrated  source  intensity  ;  the 
source  energy  is  here  assumed  to  be  1.25  MfeV  with  no  degradation. 

Changing  varietbles,  ttie  relative  exposure  due  to  an  infinite  plane 
source  is 


D 

D 


Fee' 


P(d) 


d(sin9> 

sin9 


3.2.5 


To  find  the  effect  of  source  anisotropy,  the  ratio 


F(e:  p  d'< 


d  t  sin9) 
sin6 


P(d> 


d  i  sj  n9, 
sin9 


was  calculated  and  found  to  be  0,96  . 

This  calculation  shows  that  if  iree  field  measurements  were  acne  by 

(7) 

the  method  described  by  Fexroad  and  Schmoke  using  a  source  similar  to 

<81 

that  described  by  Burson  and  Summers  ,  source  anisotropy  would  decrease 
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the  measured  relative  exposure  by  approximately  4%  .  If  the  source  were 
such  that  the  exposure  on  the  axis  were  greater  than  that  observed  normal 
to  the  axis,  source  anisotropy  would  increase  the  relative  exposure. 

3,3  Effect  of  Lack  of  Source  Reflection  in  Roof  Experiments 

The  barrier  and  geometry  factors  in  present  methods  of  calculation 
of  fallout  protection  are  calculated  usually  by  the  method  of  moments. 
These  assume,  among  other  things,  an  infinite  extent  of  material  on  that 
side  of  the  source  which  is  away  from  the  detector. 

In  reality,  fallout  would  be  in  a  thin  layer  on  and  around  a  struc¬ 
ture  with  air  above  it.  If  the  fallout  field  were  infinite  in  extent, 
or  greater  than  three  mean-free-path  lengths  (about  1500  feet)  in  radius, 
and  if  the  building  height  were  small  compared  with  a  mean-free-path 
length  in  air,  then  the  contribution  of  radiation  scattered  from  the  air 
above  would  equal  the  theoretical  contribution  due  to  the  radiation 
scattered  back  by  the  hypothetical  infinite  medium. 

If  the  source  in  a  roof  experiment  were  covered  with  an  infinite 
thickness  of  the  same  material  as  the  roof,  the  experimental  configura¬ 
tion  would  be  identical  with  the  theoretical  source  configuration  in  the 
moments  calculations, and  this  would  include  the  equivalent  of  "skyshine." 
Since  this  was  not  always  the  case  in  past  experiments,  an  estimate  must 
be  made  of  how  much  the  theory  may  have  differed  from  the  experimental 
results  on  account  of  the  experimental  source  characteristics. 

(5) 

Berger  and  Doggett  calculated  buxld_up  factors  for  finite  and 
semi-infinite  barriers  for  plane  normal  sources.  The  difference  between 
the  results  for  the  two  barriers,  for  a  1  MeV  source  on  iron,  at  one 
mean-free-path,  was  1%  .  For  an  isotropic  source,  the  difference  may 
be  expected  to  be  larger. 

An  approximate  calculation  can  be  made  as  follows:  From  the  results 
presented  in  Appendix  A  the  exposure  penetrating  into  iron  from  a  i.25 
MeV  plane  isotropic  source  is  known.  We  wish  to  subtract  the  exposure 
due  to  the  radiation  which  moves  initially  away  from  the  detector  and 
then  cresses  the  source  plane  to  be  detected. 
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SUMMARY 


A  comparison  is  made  between  experiment  and  theory  for  the  calcu¬ 
lation  of  shielding  by  a  structure  against  fallout  radiation. 

Experimental  results  for  exposure  penetrating  a  roof  slab,  and 
for  reduction  factor,  are  found  to  agree  with  moments  theory  calcu¬ 
lations,  often  to  better  than  10%  when  the  geometry  factor  L  (u),X) 
was  used. 

A  comparison  between  experiment  and  theory  may  be  inaccurate  due 
to  anisotropy  of  the  experimental  source,  to  lack  of  source  reflection 
in  root  experiments,  and  due  to  error  in  estimating  the  thickness  of 
the  roof.  The  magnitudes  of  the  errors  due  to  these  effects  are  inves¬ 
tigated,  and  found  to  be  small,  but  not  necessarily  negligible. 

Detailed  results  of  penetration  in  iron  and  concrete,  due  to  plane 
Isotropic  sources  of  various  energies  are  given,  and  a  review  of  ground 
contribution  to  fallout  radiation  penetration  through  simple  structures 
is  appended. 
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Berger  and  Raso  give  the  following  values  for  the  number  albedo 
A  from  a  1.25  MeV  isotropic  source:  on  iron  0.295  ;  on  concrete 
0.383. 

(9) 

Berger  and  Raso  calculated  the  reflected  energy  spectrum  ior  a 
1  MeV  isotropic  source  on  iron.  This  is  shown  in  Fig.  3.  On  the  basis 
of  this  result,  i‘  will  be  assumed  that  all  the  reflected  photons  due  to 
a  Co-60  source  on  iron  or  concrete  have  energy  0. 25  MeV  . 

The  current  due  to  scattered  radiation  which  crosses  the  source 
plane  will  be  A  times  the  source  current  and  is  assumed  to  be  at  0.25 

MeV  .  Assuming,  without  justification,  that  this  is  isotropic  at  the 
source  plane,  the  penetration  of  the  scattered  radiation  can  be  found 
from  the  moments  calculation  at  0.25  MeV  .  The  exposure  at  depth  Z 
psf  in  the  barrier,  due  to  reflection  above  the  source,  will  be 

A  •  D(Z,  0.25  MeV) 

D(Z,  1.25  MeV) 

This  calculation  will  be  approximate,  because  the  reflected  photon  dis¬ 
tribution  will  be  more  peaked  than  isotropic  and  because  the  energies 
of  the  reflected  photons  are  not  at  a  single  energy,  but  an  estimate  of 
the  order  of  magnitude  of  the  effect  is  obtained  from  this  calculation. 
The  calculated  ratios  of  exposure  reflected  by  source  cover  to  exposure 
in  an  infinite  medium,  are  shown  in  Table  2. 

The  above  calculation  gives  approximately  the  same  result  as  cal¬ 
culations  using  the  functions  given  by  Spencer^ \  In  his  notation: 

S<"d)  is  the  ratio  of  exposure  due  to  radiation  from  the  upper  hemisphere 
(  akyshine")  to  the  reference  exposure  3  feet  in  air  above  the  source 
plane; 

S' (X)  is  the  ratio  of  exposure  at  X  psf  to  the  exposure  at  0  psf,  due  to 
a  plane  source  in  an  infinite  medium,  emitting  radiation  Isotropically 
into  the  half-space  in  which  the  detector  is  buried; 

L(X)  is  the  rsito  of  exposure  due  to  a  plane  infinite  isotropic  source, 
at  X  psf  in  an  infinite  medium,  to  the  reference  exposure  3  feet  in  air 
above  the  source  plane. 
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A  comparison  between  the  two  methods  of  calculation  is  obtained  by 

.  A  •  D(Z,  0.25  MeV)  ,  _  ,  ,  ,  _ 

comparing  the  value  of  -  — :  ~rr  - in  the  previous  calculation 

r/h's  R«  ✓y'*  D^z>  1,25  MeV' 


with 


S(d)  S'(X) 
L(X) 


in  this  calculation,  when  d  =  0  .  Note  that  both 


calculations  refer  to  infinite  source  planes;  solid  angle  fraction  uo 


is  then  one  and  S  (X,uj)  =  1  . 

a 

The  latter  approach  was  used  by  Raso  and  Woolf 


(17) 


when  comparing 


their  Monte  Carlo  calculations  of  roof  contributions  with  the  moments 
calculations  of  Spencer^.  When  the  quantity  S(d)  S' (X)  S  (<ju,d)  was 

H 

subtracted  from  the  value  of  L(X)  L  (X,u))  ,  good  agreement  was  found 

& 

with  the  Monte  Carlo  results. 


TABLE  2 


Ratio  of  exposure  reflected  by  source  cover  to  total  exposure  for  iron 
and  concrete.  Results  from  approximate  calculation,  using  albedo  and 
penetration  curves;  and  using  curves  of  NBS  42^'. 


Material 

Z  PSF 

Calculating 

Procedure 

25 

50 

100 

150 

Iron 

0.025 

0.011 

0.003 

0.0005 

mm 

D(Z,  1.25  MeV) 

Concrete 

0.062 

0.037 

0.013 

0.004 

Concrete 

0.077 

■ 

0.031 

— 

S(d)  S*(X) 

L(X)  :  d=0 

The  results  of  the  two  mett  of  calculstion  are  given  for  concrete 
in  Table  2  and  show  reasonable  agreement  at  low  barrier  thicknesses,  when 
the  correction  is  more  significant;  at  greater  thicknesses  agreement  is 
not  as  good,  but  the  correction  is  then  not  important. 

The  calculations  show  that  lack  of  source  reflection  may  account  for 
about  3%  of  the  exposure  in  iron  and  6%  in  concrete  at  25  psf,  about  1% 
in  iron  and  4%  in  concrete  at  50  psf.  and  less  at  greater  distances.  This 
difference  between  theory  and  experiment  may  be  avoided  simply  by  having  a 
reflector  over  the  source. 
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3.4  Errors  in  Estimating  Roof  Thickness 

During  the  course  of  this  investigation,  a  question  arose  as  to  the 
error  provoked  when  relating  theoretical  to  experimental  exposures,  by 
assuming  that  a  slab  has  exactly  the  nominal  thickness  when  in  fact  it  may 
not  have.  An  extimate  of  this  cause  of  error  is  made  as  follows. 

Consider  a  roof  slab,  upon  which  is  placed  a  uniformly  distributed 
plane  isotropic  source.  The  radiation  exposure  within  the  slab  can  be 
approximately  given  by  D  =  D^e  ^  when  t  is  the  slab  thickness  and 
li  is  the  effective  broad  beam  attenuation  coefficient  for  the  slab.  The 
relation  between  errors  in  D  and  t  is  given  by  — ^  fit  .  3.4.1 

The  effective  broad  beam  attenuation  coefficient  for  a  plane  iso¬ 
tropic  Co-60  source  on  steel  is  found  from  Fig.  A2  of  Appendix  A  to  be 
0.027  psf  1  ,  or  1.1  in.  1  Assuming  this  is  known  with  zero  error 
— =  -1.1  6t  where  t  is  in  inches.  3.4.2 

The  difference  between  real  and  nominal  thicknesses  is  demonstrated 
by  observations  made  during  an  unpublished  experiment  at  the  University 
of  Illinois.  Sheets  of  glass,  graphite,  aluminum,  and  lucite,  each 

4.5  sq.  ft.  in  area  and  nominally  l"  thick,  were  measured  on  a  6"  sq. 
grid.  The  mean  dimensions  found  are  shown  in  Table  3. 

TABLE  3 

Measured  thickness  of  nominal  i ''  muterial.  The  mean  of  18  micro¬ 
meter  readings  is  given 


Materia] 

Thickness  Cinches) 

A 1 um i num 

0.980 

Class 

1.005 

Luci*e 

0.967 

Graphite 

1.002 

The  American  and  Iron  Steel  l 
for  structural  plates  to  be  minus 


(12‘* 

nstitute  gives  the  rolling  tolerances 
0%  plus  7%  .  If  we  assume  that  steel 
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used  in  an  experiment  Is 
then  6t  =0.03  x  t  and 

-jp~  =  1.1  x  .03  x  t 
=  0,03  x  t 


3%  thicker  than  the  assumed  nominal  value, 
the  observed  exposure  will  be  decreased  by 


3.4.3 


The  observed  exposure  will  then  decrease  by  about  3%  for  each  inch 
thickness  of  steel. 


4  COMPARISON  OF  THEORY  AND  EXPERIMENT 

(3) 

4. 1  The  Nuclear  Defense  Laboratory  Blockhouse  Experiment 

In  this  experiment  sources  of  Co-60  and  Cs-137  were  placed  on 
the  roof  of  a  blockhouse  so  as  to  simulate  a  plane  source;  detectors  were 
placed  within  the  blockhouse  and  exposures  ure  measured.  In  a  separate 
experiment,  the  same  sources  were  placed  »  en  field  and  the  exposure 

three  feet  above  the  effectively  infinite  ,  source  measured.  Dividing 
the  former  exposure  by  the  latter  gave  the  reduction  factor. 

The  blockhouse  had  internal  dimensions  12  ft  by  12  ft  by  8  ft  high. 
The  walls  were  of  concrete,  4  inches  thick.  The  roof,  which  was  supported 
by  the  walls  and  by  a  joist  at  the  blockhouse  center  line,  consisted  of 
1/2  inch  plywood,  and  various  thicknesses  cl  steel  or  concrete  above  this. 
The  source  rested  on  a  1/2  inch  thickness  of  lucite  on  the  roof. 

The  detectors  were  Victoreen  Model  208  lOmr  dosimeters.  These  were 
calibrated  to  ±3%  against  a  Victoreen  Model  130  r-meter  which  itself 
had  been  calibrated  to  ±2%  at  the  National  Bureau  of  Standards.  Twenty- 
four  detectors  were  placed  within  *he  blockhouse  at  various  positions. 

The  sources  used  were  0.34  curios  of  Co-€0  and  1.35  curies  of 
Cs-137.  The  cobalt  was  metal  in  a  stainless  steel  capsule;  the  cesium 
was  CsCl  also  in  a  stainless  steel  capsule.  The  sources  were  calibrated 
by  holding  them  a  known  distance  from  NBS-calibrated  ion  chambers.  The 
calibration  of  the  ion  chambers  was  accurate  to  ±2%  ,  and  it  is  expected 
that  the  calibration  of  the  source  relative  to  the  ion  chambers  was 
accurate  to  about  2%  ,  so  the  source  calibration  was  probably  accurate  to 


(0.02)2  +  (0.02) 2 


3% 
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Each  source  was  placed  on  the  roof  in  a  number  of  locations  and  the 
exposure  rates  thus  observed  were  summed  so  as  to  be  equivalent  to  a 
uniform  source  distribution  on  the  roof. 

The  experimental  exposures  at  various  heights  on  the  blockhouse 
center  lines  were  compared  with  theoretical  values  by  the  NDL  workers, 
in  the  following  way: 

(a)  The  exposures  observed  with  a  steel  roof  were  multiplied  by  a 
factor  of  1.11  to  make  them  equivalent  to  those  which  would  be  observed 
with  a  concrete  roof.  The  factor  1.11  was  used  as  the  average  ratio 

of  ordinates  of  the  curves  of  exposure  vs  effective  mass  in  concrete 
and  iron  for  a  1  Me»  plane  normal  source  (Fig.  22.1  of  reference  1). 

(b)  The  reduction  factor  was  calculated  by  dividing  this  exposure 

rate  by  the  exposure  rate  D  ,  observed  three  feet  above  an  infinite 

o 

contaminated  field.  The  values  of  D  were  observed  in  previous  work 

(7)  ° 

at  the  same  site. 

(c)  The  theoretical  reduction  factor  was  taken  from  Fig.  28.19  of 

NBS  42^\  which  gives  values  of  L(X)  L  (X,uj)  calculated  for  water  and 

a 

a  1.12  hr  fission  source. 

Comparison  of  those  experimental  results  with  theory  was  repeated 
in  this  investigation,  making  use  ot  the  newly  compu’ed  penetrations  in 
steel  (see  Appendix  A'  and  using  a  procedure  as  follows: 

(a)  Roof  thicknesses  were  expressed  as  actual  pounds  per  square  foot 
Z.  not  as  equivalent  mass  thickness  X  i see  Appendix  A). 

tb'  The  exper. mental  exposure  rates  for  steel  roofs  given  in  Tableo 

2 

3.1  and  3.2  of  the  '*DL  report  in  roentgen/hr  per  curie/ft  were  divided 
by  1.11,  so  as  to  eliminate  the  correction’'  previously  made.  The  sources 
had  been  calibrated  using  specific  gamma  ray  constants  of  14.3  roentgen/hr 
at  one  foot  from  a  curie  of  Co-63  and  4.2  roentgen/hr  at  one  foot 
from  a  curie  of  Cs-13  7,  '’he  values  corresponding  to  86.9  roentgen 
per  erg/gw  are  14.0  roenigen/hr  at  one  foot  from  a  curie  of  Co-60  and 
3.4  roentgen/hr  at  one  foot  from  Cs-137.  The  values  of  exposure  rate 
per  curie,  used  in  the  theory  and  injthe  reduction  of  experimental  data, 
were  made  equal. 
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(c)  The  exposures  obtained  from  the  moments  calculation  for  1.25 

2 

MeV  were  multiplied  by  5284  to  convert  from  MeV/gm  per  photon/cm  to 

2 

roentgen/hr  per  curie/ft  .  This  constant  is  obtained  as  follows: 

1  MeV  =  1.602  x  10"6  erg 
1  roentgen  =  86. 9  erg/gm 

1  curie  of  cobalt  =  7.4  x  10^  photons/sec 

2  2 
1  ft  =  929.03  cm 

2 

Thus,  1  MeV/gm  sec  per  photon/cra  sec 


1.602  x  10 
86.9 


-6 


1U 


7„ 4  x  10 
*  929.03 

=  5284  roentgen/hr  per  curie/ft' 


x  3600 


Since  cesium-137  emits  0.84  photons  per  disintegration  whereas 

cobalt-60  emits  two  photons  per  disintegration,  the  conversion  factor  for 

2 

cesium-137  will  be  2219  roentgen/hr  per  curie/ ft  . 

(d)  The  geometry  factors  L  (qa,X)  and  L  (uj,X)  were  obtained  from 

a  (1 )  C 

Figs.  B27 ,  B28,  B29,  and  B30  of  NBS  42  ,  which  are  for  Co-60  and  Cs-137 

in  concrete.  Geometry  factors  have  rot  been  calculated  for  steel,  because 
it  is  expected  that  these  factors  will  be  much  less  dependent  on  the 
difference  between  mass  attenuation  coefficients  for  steel  and  concrete 
than  the  barrier  factor. 


(e)  The  theoretical  exposure  rates  were  computed  from  D(Z)L  <a),X'' 

a 

and  D(Z)L  (tu,X)  ,  where  D(Z'  was  obtained  as  shown  in  (c)  above, 
c 

These  are  compared  with  the  NDl  observations  in  Table  4. 

(f)  The  reduction  factor  was  ca.culated  by  dividing  the  exposure 
rate  obtained  in  te)  above  by  the  calculated  exposure  in  the  air  three 
feet  above  the  ground  as  shown  in  Appendix  A.  These  calculations  are 
compared  with  the  NDL  olservations  in  Table  5. 


Tables  4  and  5  and  Figures  4  though  7  show  the  calculated  and  observed 
exposure  rates  and  reduction  factors.  The  curves  are  not  plotted  beyond 
45  paf  because  the  observations  at  50  psi  were  for  a  concrete  roof,  whereas 
the  other  observations  and  calculations  are  for  ste?l. 

There  are  several  possible  causes  of  error  in  comparing  theory  to 

(3) 

experiment;  since  these  were  not  all  evaulated  in  the  report  an  order 
of  magnitude  estimate  of  these  errors  is  made  here. 
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TABLE  4 


Ratio  of  Theoretical  Exposure  Rate 
2 

(mr/hr  per  curie/ft  )  to  Observed  Exposure  Rate 


Source 

Detector 

(11 

Barrier  Thickness  PSF 

Theoretical 
Geometry 
Factor  Used 

Location 

18. 0  + 

33. 7* 

44.  9+ 

EES ■ 

Co- 60 

C  6 

0.73 

1.08 

1.08 

1.09 

0.98 

L  (uu,X) 
a 

C  3 

0.41 

0.97 

1.03 

1.08 

0.97 

C  0 

0.24 

0.74 

0.83 

0.92 

0.83 

C-l 

0.  20 

0.79 

0.86 

0.85 

0.78 

Co- 60 

C  6 

0.73 

1.17 

1.26 

1.29 

1.22 

Lc(-jj,X) 

C  3 

0.41 

0.95 

1.07 

1.18 

1.17 

C  0 

0.24 

0.72 

0.83 

0.JH 

0.94 

C-l 

0.20 

0.  7S 

0.84 

0.91 

0.90 

i 

Cs- 137 

C  6 

0.73 

1.06 

1.12 

1.26 

1.  13 

L  (x .  X) 
a 

C  3 

0.41 

1.06 

1.02 

1.06 

1.03 

C  0 

0.24 

0.97 

0.91 

0.92 

0.90 

C-l 

0.20 

0.98 

0.87 

0.91 

0.94 

Cs-137 

C  6 

0.73 

1.26 

1.37 

1.55 

1.45 

L  (x.Xl 
c 

C  3 

0.41 

0.96 

1.25 

1.35 

1.43 

C  0 

0.25 

0.81 

1.03 

1.11 

1.22 

■  — . 

C-l 

0.20 

0.78 

0.97 

1.06 

1.18 

*  Steel  was  used  in  the  experiment  for  these  thicknesses. 

*  Concrete  wss  used  in  the  experiment  for  this  thickness  only. 
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TABLE  5 


Ratio  of  Theoretical  Reduction  Factor 
to  Observed  Reduction  Factor 


Source 

Detector 

Location 

(D 

Barrier  Thickness  PSF 

Theoretical 

18.  0+ 

33.  7+ 

44.  9+ 

50.3* 

Geometry 
Factor  Used 

Co-60 

C  6 

0.73 

1.00 

1.00 

1.00 

1.02 

L  (u),X) 

C  3 

0.41 

0.  82 

0.95 

1.  00 

1.03 

C  0 

0.24 

0.69 

0.77 

0.84 

0.85 

C-l 

0.20 

0.71 

0.74 

0.-81 

0.78 

C0-60 

C  6 

0.73 

1.09 

1.17 

1.20 

1.26 

L  (u),X) 

c 

C  3 

0.41 

0.80 

0.98 

1.09 

1.  22 

C  C 

0.24 

0.67 

0.77 

0.91 

0.96 

C-l 

0.20 

0.71 

0.  74 

0.85 

0.91 

Cs-137 

C  6 

0.73 

0.94 

1.00 

1.  13 

1.10 

La(u),X) 

C  3 

0.41 

0.94 

0.91 

0.97 

0.99 

C  0 

0.24 

0.87 

0.81 

0.82 

0.87 

C-l 

0.20 

0.  86 

0.79 

0.84 

0.94 

Cs-137 

C  6 

0.73 

1.12 

1.  2? 

1.39 

1.42 

Lc  (tu,X) 

C  3 

0.41 

0.86 

1.11 

1.20 

1.42 

C  0 

0.  24 

0.71 

0.91 

1.00 

1.17 

C-l 

o.?o 

0.  70 

_ 

0.86 

0.95 

1.17 

Steel  was  used  in  the  experiment  for  these  thicknesses, 

* 

Concrete  was  used  in  the  experiment  for  this  thickness  only. 
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(a)  The  error  ia  the  calculated  values  was  of  the  order  of  2%  due 

to  difficulties  in  reading  accurately  the  values  of  L  (u),X)  and  L  (d,X) 

( 1 '  3  C 
from  the  graphs  in  N?S  42  ' ; 

(b)  Uncertainty  in  absolute  calibration  of  the  source  strength  is 
expected  to  have  been  approximately  37.  (see  text); 

(7) 

vC)  Because  of  anisotropy  of  the  source  in  free  field  measurements  , 

2 

D  for  Co-60,  when  expressed  as  exposure  rate  per  curie/ft  ,  may  possibly 
o 

have  been  in  error  by  approximately  4%  (see  Section  3.2); 

(d)  If  the  Co-60  source  were  entirely  unreflected  then  the  theoretical 
exposure  rate  would  have  been  higher  than  the  experimental  exposure  rate 

by  approximately  2%  at  25  psf  and  1%  at  50  psf  in  steel  (see  Section  3.3). 
Since  in  fact  the  finite  size  of  the  source  would  have  caused  some  radia¬ 
tion  to  be  reflected  back  into  the  roof  this  effect  may  have  been  small, 
for  the  steel  roofs  tested.  For  concrete,  the  effect  would  have  been 
larger,  but  still  small; 

(e)  The  nominal  thicknesses  of  the  steel  roofs  appear  to  have  been 
used  for  the  calculations  in  the  NDL  report.  II  is  possible  that  the 
thicknesses  were  in  fact  greater  than  nominal.  If  they  were  3%  thicker 
the  exposures  would  have  been  cvercalculated  by  approximately  1%  at  18  psf, 
2%  at  34  psf  and  3%  at  45  psf  for  the  Co-60  source  (see  Section  3.4); 

(f)  Because  of  degradation  of  the  source  spectrum,  the  measured 

2 

values  of  exposure  rate  per  curie/ft  through  the  barrier  would  have  been 
lower  than  the  theoretical  values.  ,T,he  magnitude  of  this  effect  is,  for 
the  present,  unknown  (see  footnote  to  page  6'. 

Ihe  effect  listed  in  fcj  above  would  have  reduced  or  increased  the 
free  field  exposure  rate,  depending  upon  the  geometry  of  the  encapsulated 
source,  whereas  those  listed  in  (d),  (e),  and  (f)  above  would  have  decreased 
the  exposure  rate  observed  in  the  blockhouse.  This  discussion  of  possible 
errors  is  based  on  general  considerations  and  no  judgement  can  be  given 
here  as  to  whether  ~hese  causes  of  difference  between  theory  and  experi¬ 
ment  were  actually  present. 


19 


Agreement  between  calculated  and  observed  reduction  factors  appears 

to  be  better  for  larger  solid  angles  and  greater  barrier  thickness,  and 

better  agreement  is  found  between  theory  and  experiment  when  the  geometry 

factor  L  (cd , X)  is  used  than  when  the  factor  L  (oj,X)  is  used, 
a  c 

The  geometry  factor  L  (u),X)  is  calculated  for  an  infinite  plane 

source  in  an  infinite  medium  and  a  limited  solid  angle  of  detector  response 

the  factor  L  (u),X)  is  calculated  for  a  circular  plane  source  in  an 
c 

infinite  medium  and  a  detector  response  over  4tt  steradiams.  Since  in 
fact  the  barrier  has  a  finite  thickness,  and  the  detector  is  some  dis¬ 
tance  from  the  barrier  with  only  air  in  its  close  vicinity,  the  factor 
L  (u>,X)  may  indeed  be  expected  to  give  more  accurate  results. 

Of  the  16  values  of  reduction  factor  calculated  for  Co-60  (using 

L  (u),X)  ): 

8 

7  fall  within  ±  10%  of  experiment; 

11  fall  within  ±  20%  of  experiment; 

15  fall  within  ±  30%  of  experiment. 

For  Cs-137  (using  L  (u),X)  ): 

8 

8  fall  within  ±  10%  of  experiment; 

15  fall  within  ±  20%  of  experiment. 

(2)* 

4.2  The  Protective  Structures  Development  Center  Idealized  Roof 

In  this  experiment,  a  circulating  source  was  placed  on  the  ground 
with  concrete  slabs  above.  Detectors  were  placed  at  various  heights  above 
the  slabs. 

The  reduction  factors  were  calculated  by  dividing  the  observed  expo¬ 
sure  rates  by  the  exposure  rates  previously  measured  three  feet  above 
an  infinite  plane  of  contamination.  These  experimental  reduction  factors 
were  compared  with  the  theoretical  values  of  L(X)L  (X,ua)  and  L(X)L  (X,u>) 

A  C 

and  this  comparison  is  repeated  here  by  taking  the  values  given  in  that 

♦  (4) 
report. 


* 


At  the  time  of  this  writing,  preliminary  results  only  are  available. 
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The  L(X)L  (X,uj)  curves  were  closer  to  the  observed  results  than 
a 

the  L(X)L  (X ,uj)  curves,  at  all  values  of  X  .  The  differences  are 
c 

listed  in  Table  6, 


TABLE  6 


Comparison  of  Theory  and  Experiment  for  the  PSDC  Idealized  Roof. 
Ratio  of  Calculated  Reduction  Factor  to  Observed  Reduction  Factor 


Barrier 

PSF 

U) 

Theoretical  Value 

L(X)L  Cut)  i  X) 
a 

L(X)L  (uu.X) 

C 

0  | 

0.02 

■ . 

1.00 

0.95 

0.  10 

1.00 

0.94 

0.60 

1.00 

0.94 

48 

0.02 

1.12 

1.19 

0.10 

1.15 

1.25 

0.60 

1.17 

1.43 

96 

0.  02 

1.07 

1.40 

0.  10 

1.13 

1.47 

0.60 

1.09 

1.36 

144 

0.02 

1.00 

1.50 

0.  10 

1.00 

1.55 

0.  60 

1.00 

1.35 

In  this  experiment,  the  unprotected  exposure  was  observed  with  detec¬ 
tors  above  the  circulating  source.  The  protected  exposure  was  obtained 
by  inserting  concrete  slabs  between  source  and  detectors.  In  this 
arrangement,  error  duo  to  lack  ol  source  reflection  was  eliminated  since 
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the  ground  served  as  a  reflector.  The  source  was  a  cylinder  which  cir¬ 
culated  so  that  the  circular  suiface  faced  the  dosimeters  for  a  large 
portion  of  the  travel.  Since  source  anisotropy  is  usually  most  marked 
beyond  the  ends  of  the  source,  it  is  expected  that  this  arrangement 
gave  less  error  due  to  source  anisotropy  than  the  previously  discussed 
experiment  where  a  single  source  was  placed  on  the  ground  or  on  the  roof. 

Because  of  the  two  above-mentioned  reasons,  and  because  the  experi¬ 
ments  w.  re  compared  with  calculations  for  Co-60  incident  on  concrete 
slabs,  th^-.e  experimental  results  were  compared  with  existing  theo¬ 
retical  calculations,  no  further  calculations  being  carried  out. 

5  CONCLUSIONS 

Various  possible  causes  cf  error  in  comparing  experiment  with  theory 
were  onsidered  in  this  investigation.  The  conclusions  may  be  summarized 
as  follows: 

(1)  Source  anisotropy  in  free  field  measurements  using  a  stationary 
source  may  cause  an  error  on  the  order  of  4%  when  normalizing  the 
measured  exposure  to  the  calibrated  source  intensity.  It  must  be  empha¬ 
sized  that  this  is  an  approximate  result  only.  If  more  accuracy  in 
estimating  this  cause  of  error  is  required,  more  detailed  calculations 
must  be  made  for  any  particular  experiment. 

(2)  The  lack  of  reflector  above  a  source  in  a  roof  experiment  will 
cause  larger  error  for  a  concrete  than  for  a  steel  roof;  this  error  for 
concrete  will  be  of  the  order  of  6%  at  25  psf,  3%  at  50  psf,  and 
decreases  to  less  than  1%  above  100  psf.  The  source  should  have  a 
reflector  so  as  to  include  the  equivslent  of  "skyshine"  which  is  implicit 
in  the  moments  theory  calculations. 

(3)  Errors  in  the  assumed  value  of  thickness  of  a  steel  roof  will 
cause  an  error  in  the  calculated  exposure  of  about  1%  per  inch  thick¬ 
ness,  for  an  error  of  1%  in  assumed  roof  thickness. 

(4)  Comparison  of  the  results  of  experiments  with  theoretic;  1  values 
often  show  agreement  to  within  10%  when  the  geometry  factor  L^fui.X)  is 
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used.  There  is  a  trend  for  the  theory  to  predict  lower  values  of  expo¬ 
sure  and  reduction  factor  as  solid  angle  decreases. 

(5)  It  is  recommended  that  the  geometry  factor  La(uu,X)  be  used 
for  roof  penetration  calculations  in  simple  structures,  and  in  that  case 
an  accuracy  of  greater  than  10%  in  predicted  reduction  factor  may  often 
be  expected. 
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Enlarged  view 
of  source 


SCHEMATIC  VIEW  OF  FREE  FIELD  EXPERIMENT  AT 
THE  NUCLEAR  DEFENSE  LABORATORY  (REFERENCE  7) 


FIG.  2 


2 


FIG.  3  ENERGY  SPECTRUM  OF  I  MeV  PHOTONS 
REFLECTED  FROM  IRON 


ISOTROPIC  INCIDENCE,  I  MeV  SOURCE 
FROM  TABLE  A  19,  BERGER  S  RASO 

NBS  5982,  I960 


FIG. 4  COMPARISON  OF  EXPERIMENTAL 
AND  THEORETICAL  EXPOSURE 
RATES 


EXPERIMENTAL  —  SCHMOKE  AND  REXROAD  1962 

o  Detector  Position  C6 

•  C3 

a  CO 

•  C-l 

Theoretical  -  La(w,X)D(X) 

- Lc  (w,X)  D  (X) 


Source  Co-60 
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FIG. 5  COMPARISON  OF  EXPERIMENTAL 
AND  THEORETICAL  EXPOSURE 
RATES 


EXPERIMENTAL  -  SCHMOKE  AND  REXROAD  1962 

o  Detector  Position  06 

•  C3 

*  CO 

»  C-l 

L„  (asX)D(X) 

Lc  (u>,X)  D  (X) 

Source  Cs-137 


Theoretical 


o 


C6 


REDUCTION  FACTOR 
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BARRIER  PS  F 


FIG. 7  COMPARISON  OF  EXPERIMENTAL 
AND  THEORETICAL  REDUCTION 
FACTORS 


EXPERIMENTAL  -  SCHMOKE  AND  REXROAD  1962 

o  Detector  Position  C6 
•  C3 

a  CO 

■  C-l 

Theoretical  - -  La  (<u,X)  D  (X)/D0 

— - Lc  («,X)  D  (X)  /D0 

Source  Cs-137 
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APPENDIX  A 

CALCULATIONS  OF  EXPOSURE  AND  REDT.'CTION  FACTOR  DUE  TO 
PLANE  ISOTROPIC  SOURCES  IN  IRON  AND  CONCRETE 


Exposure  distributions  were  calculated  in  infinite  media  of  concrete, 
iron  and  air  for  plane  isotropic  sources,  by  the  method  of  moments.  Use 
was  made  of  the  code  available  at  the  National  Bureau  of  Standards  for 
this.  Source  energies  were  1.33,  1.25,  1.17,  0.6616,  0.25,  0.10,  and 
0.05  MeV.  The  material  cross  sections  were  taken  from  those  published 
by  Grodstein^*^  and  McGinnies^^. 

Table  A1  gives  the  exposure  in  iron  and  in  concrete  at  a  number  of 
depths  for  the  various  source  energies.  These  values  are  plotted  in 
Figures  A1  to  A4. 


Table  A2  shows  the  following  quantities  for  each  of  five  source 

energies;  the  exposure  in  air  three  feet  from  the  source  plane  in  MeV/gm 
2 

per  photon/cm  ,  and  the  ratio  of  scattered  to  total  exposure  at  the 
detector. 

2  2 

Conversion  from  MeV/gra  per  photon/cm  to  roentgen/hr  per  curie/ft 
is  shown  in  Section  4.1.  Reduction  factors  may  be  obtained  by  dividing 
the  exposures  in  concrete  or  iron  by  the  exposure  in  air  three  feet 
from  the  source  plane  at  the  relevant  source  energy. 


It  should  be  noted  that  the  exposures  are  given  as  functions  of  the 

mass  thickness  Z  .  This  is  not  to  be  confused  with  the  equivalent 

(1  ^ 

mass  thickness  X  of  NBS  42 
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TABLE  A1 


Exposures  in  MeV/gm  per  photon/cm  for  Iron  and  Concrete 
at  Different  Source  Energies  and  Depths 


E 

o 

MeV 

p  z 
o 

Iron 

Concrete 

z 

psf 

Exposure 

z 

psf 

Exposure 

1.33 

0.5 

19.9 

0.0185 

18.6 

0.0199 

1.0 

39.7 

0.00932 

37.1 

0.0103 

1.5 

59.6 

0.00520 

55.  7 

0.00588 

2.0 

79.5 

0.00302 

74.2 

0. 00347 

2.5 

99.3 

0.00179 

92.8 

0.00209 

3.0 

119 

0.00107 

111 

0.00127 

4„  0 

159 

0. 000395 

148 

0.000473 

8.0 

318 

0. 00000779 

297 

0.00000769 

1.25 

0.5 

19.3 

0.0178 

18.0 

0.0193 

1.0 

38.5 

0.00901 

35,9 

0.0101 

1.5 

57.8 

0. 00504 

53.9 

0.00575 

2.0 

77.0 

0. 00293 

71.9 

0  00340 

2.5 

96.3 

0.00174 

89.4 

0.00205 

3.0 

116 

0.00105 

108 

0.00125 

4.0 

154 

0. 000387 

144 

0. 000469 

8.0 

308 

0.00000771 

288 

0. 00000973 

1.17 

0.5 

18.6 

0.0171 

17.4 

0.0186 

1.0 

37.2 

0.00867 

34.8 

0.00973 

55.  9 

0.00486 

52.2 

0.00559 

74.5 

0.00284 

69.5 

0. 00332 

93.1 

0.00169 

86.9 

0.00201 

3.0 

112 

0.00102 

104 

0.00123 

4.0 

149 

0.000376 

139 

0. 000463 

8.0 

298 

0.00000759 

278 

0. 00000977 

.6616 

0.5 

14.  1 

0.0113 

13.3 

0.0130 

1.0 

28.2 

0.00583 

26.5 

0. 00709 

1.5 

42.  2 

0.00332 

39.8 

0.00421 

56.3 

0.00197 

53.0 

0.00258 

70.4 

0.00119 

66.3 

0.00160 

84.4 

0. 000723 

79.5 

0.00100 

4.0 

113 

0.000274 

106 

0. 000395 

8.0 

225 

0.00000597 

212 

0.00000961 

34 


TABLE  A1  (Continued) 


0.25 


0.10 


0.05 


Iron 

Concrete 

z 

psf 

Exposure 

......  .  .. 

z 

psf 

Exposure 

8.66 

0.  00400 

8.93 

0.00584 

17.3 

0.  00200 

17.9 

0. 00335 

26.0 

0.00111 

26.8 

0.00206 

34.  7 

0.000636 

35.7 

0.00130 

43.3 

0.000374 

44.6 

0. 000822 

52.0 

0.000223 

53.6 

0. 000525 

69.3 

0.0000806 

71.4 

0.000215 

139 

0.  00000155 

143 

0. 00000588 

3.00 

0.000901 

6.03 

0.00193 

5.99 

0.000389 

12.1 

0.00105 

8.99 

0.000189 

18.1 

0.000613 

12.0 

0.  0000978 

24.1 

0. 000366 

15.0 

0.  0000522 

30.1 

0. 000221 

18.0 

0.  0000286 

36.2 

0. 000135 

24.0 

0. 00000893 

48.2 

0.0000507 

47.9 

0.000000110 

96.4 

0. 00000108 

0.569 

0. 000577 

2.92 

0. 000972 

l.H 

0.000232 

5.84 

0.000448 

1.71 

0.000107 

8.76 

0. 000228 

2.28 

0.0000530 

11.7 

0.000121 

2.85 

0.  0000273 

14.6 

0. 0000665 

3.42 

0.0000144 

17.5 

0.0000371 

4.55 

0.00000425 

23.4 

0.0000120 

9.11 

0. 0000000443 

46.7 

0.000000163 
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TABLE  A2 

Exposures  in  Air  Three  Feet  from  a  Plane  Isotropic  Source 


E 

o 

MeV 

Exposure 

2 

MeV/gm  per  photon/cra 

_  ^  Scattered  Exposure 

Ratio  — - - - 

Total  Exposure 

1.33 

0.0957 

0.172 

1.25 

0.0911 

0.177 

1.17 

0. 0865 

0.181 

0.6616 

0.0530 

0.231 

0.05 

0.00624 

0.503 

From  the  above  table,  it  can  be  calculated  that: 

For  Co-60,  exposure  rate  three  feet  above  source  plane  a  482  roentgen/hr 
2 

per  curie/ft  ; 

For  Cs-137,  exposure  rate  three  feet  above  source  plane  =  118  roentgen/hr 
per  curie/ft^  . 


The  assumed  density  of  air 


0.001205  g/cc  . 


37 


U 


=“  =  —  - 


!====  = 


r=r-rr^  =  =:  =  ==-==  =  : 


ii^llliiigiis 


FIG.  A2  EXPOSURE  (  MeV/gm  per  photon /cm  ) 

vs  DEPTH 


■■IHMI 


E0=  1.25  MeV  8  .6616  MeV 

- Iron 

-  Concrete 

Source- Plane  Isotropic 


■■■■■■■■■■■■■■■■■■■ 


S  =  S 


_  1 

i  *  10' 


llli 


illiS 


==,==== 


===•=== 


============== 


t^iil 


rSSssii 


iSS 


Ii==isillslll2=gllilil= 


•  rr  —  -  - - r  — 


=  ==  -  •  = 


lisiiiii 

:=5S55S| 


ill! 


liiS 


IS 


liiililiiHiiiaiiiiHSlHiKBil 

■Kiii&iiiHifiiili&l 


iHk  =.«**»  w  r-  i 


?  =  S:  =  =  |  =  i  =  i  =  SE=  =  sii=25 

E=ii=====i=ii=ii=i=I 
llllllllllil 


22;: 

-=s===========,=========i^====-— ™  »iss‘?ai'^ssssii? 

=======§ =11=1  i|i  11  ||!§i  slip  sjjjsilfl  iisssi‘3§.“s =s=5i= I 


ssEs:sisss8ssiiiiiiisisjjs|sssissSsss:2ssissiis|iss|Ksj|^p^»^s 

■■gggggggggggggggggg  gggggSSSSSS55SSSSSg5SS5SSSS5S55S8««  8SS85 


liHiiiiiiiiliiiiiliipiii! 

-  .  - _  I^Siiisi  j4EsSS**MMM^ 

; is|siis|iiii3lliiiif 

—  ==?*§i*s£SsSssiSSSfi£i 

illlfipiiiiiriisa 

IflllMllilllilil 

immm 


I  *  to 


-  =3*=£siii3=2&!sl=§Misi 


L = iiiiiii^ii^liiiiiill 

|||[j|5ss5axig£i:5|»|5||l 

ittaifeBP 


i « io‘ 


IliiiiiglliSliigSaii^j 

|ss 5  ESE===iHMMi 


iiiii iiii!iiii| 


■Kill 

iSlSiiSii, 


=  =  ■=>■- jS  I 

SISSSI 


-  »  .  ,'b;^  ■>  2, 

-f  .ii.  :  „  %-mAid'Ssm 

;  =  3s  =  ss= sss&saM&Mtgtem 


I=rsa5s-ssss*si 

-  ..  _  ••  _  •■•: .  -  sm.  Sr  a* 

;  =  *££■  =iS’il£«SSa 
^si|s  liaiCCS ill 


.&  dL  ci  3k;  Ss 

"**  ‘  ~  *  ££  i*  ifi  » j*i  v-'  5l  is  £  &  ft-ifc 

^esfiBSti 


A^itbtu  Jk JSt 

a  S  a  4’  it4 ')*  ai  i«  a  s*  it sa  jbes  a>  a  aLss-ss  as  at  a  «ts 


rtl  *6*  ^«S29(f  =*  -.3#  33t  in  -*  ii  *  s»  ^  ii  *  »  ii  a»a  i- SS-t,  i»  I 


**ssss  s*53B5v isS9Ek«aii«  ait  a  »»«**«*)  autif  £•••«! 

" 'ssasaasLass  ?a*s-3323 

illltlftfltfllf  li 


zZ9t  Js  4k-  zn  ^4  4* 


2;  S=j  .^.  rr  «*  ,z&  *  a  C’4  •  a  =j  a  n*£:s:fcs!fc'aii*«*3*a§j*:a*ffl*:»isa«^siiXfBsfci!5i«g£iRjf!i| 


I_i\;3si-iis£5*'=5ES3?asa#'aBssi|iSs  =  si;s3«3iE3s®#3^323li:5s|i4a3 

r  =  k=  ai:£«£siSSSa»#*SfefcS^&**S«;#a#*#*«»£#*S*ai**fc*****«*iS**!H?l«Sfc 


jagaiaassj jfff: 


iillisHif 


40 


APPENDIX  B 

A  REVIEW  OF  GROUND  CONTRIBUTION  TO  FALLOUT  RADIATION 
PENETRATION  THROUGH  SIMPLE  STRUCTURES 


B.  I  Introduction 


There  has  been  a  flurry  of  experimental  efforts  in  the  past  few 
years  to  check  the  validity  of  the  functions  used  in  computing  the  ground 
contamination  contribution  to  fallout  radiation  penetration  into  simple 
structures,  under  infinite  field  conditions.  In  many  cases,  it  is  dif¬ 
ficult  to  design  experiments  which  will  isolate  the  function  to  be 
studied;  and  therefore  the  results  obtained  do  not  always  lead  to  unam¬ 
biguous  conclusions.  It  is  believed  appropriate  at  this  time  to  review 
some  of  the  basic  technology  and  the  conceptual  basis  therefor.  This 
will  provide  an  insight  into  the  validity  of  previous  experimental  pro¬ 
grams  and  will  be  a  basis  for  recommendations  as  to  the  directions  for 
fruitful  research  in  the  near  future. 


Let  us  consider  the  basic  formula  presently  used  to  calculate  the 
ground  contribution  to  the  attenuation  function  of  a  simple,  rectangular- 
one-story  structure  under  inf  in;: v  fallout  field  condi* :  ions: 


R  r 


B 

w 


Cd  Ed  ° 


S  > 

w 


s. 


(1) 


The  notation  used  is  roughly  similar  to  that  in  the  OCD  Engineering  Manual 

* 

and  should  be  readily  identified  bv  anyone  familiar  with  the  Manual. 

We  will  attempt  »o  explain  this  formula  in  ter*.*  of  ‘'operational'* 
concepts,  that  is,  In  terms  of  exper teen's  (perhaps  Idealized)  devised 
to  measure  them.  Concepts  which  cannot  be  explained  in  such  a  manner 


This  is  essentially  the  sameas  the  formula  lor  C  given  in  Example 
4-3  of  the  Engineering  Manual  .  Our  terms  0^  ind  arc  actually 

split  into  two  terms  in  the  Manual.  (*ur  term  Erf  does  not  appear  in 
the  Manual  formula,  but  this  is  explained  below.  The  remaining  corres¬ 
pondence  of  terms  is  obvious. 


(13) 
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cannot  be  determined  experimentally.  The  various  terms  in  the  formula 
can  be  related  to  ratios  of  exposure  rate  measurements  carried  out  foi 
various  structural  configurations.  These  are  depicted  in  Figure  B-l 
and  are  discussed  below. 


B.2  Explanation  of  Situations  Depicted 

Before  discussing  the  interrelationships  of  these  configurations, 

or  "situations,"  presented,  it  is  desirable  to  describe  them  and  to 

indicate  the  present  theoretical  (mathematical)  approach  to  finding  the 

exposure  rates  in  the  designated  detector  locations.  The  formulations 

(2) 

are  taken  or  inferred  from  Eisenhauer's  NBS  Monograph  76  ,  with  further 

reference  as  necessary  to  Spencer's  NBS  Monograph  42^\  (These  will 
hereinafter  be  called  NBS-76  and  -42,  respectively. ) 


in  terms  of  the  ^-function  defined  in  NBS-42.  This  formula  is  exact, 
and  in  fact  the  ^-function  is  normalized  so  as  to  make  equal  to 
unity. 


Situation  2  provides  an  exposure  rate  behind  two  walls,  very  close 
together.  The  wails  are  of  indefinite  height  and  length,  because  the 
detector  is  so  close  to  the  walls  that  radiation  can  be  assumed  to  be 
entering  the  detector  from  all  directions  (4tt  steradians).  The  wall  mass 
thickness  is  the  same  as  that  of  the  given  blockhouse  structure.  The 
exposure  rate  in  this  situation  is  based  on  the  NBS-42  function  for 
penetration  through  a  single  wall  (Eq.  27,7): 
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This  function  must  be  multiplied  by  2,  to  account  for  contribution  through 
two  walls  rather  than  one.  “"he  result  is  then  slightly  too  high  for 
small  values  of  the  mass  thickness,  X,  because  the  function  W  always 
includes  the  effect  of  an  infinitely  thick  backscatterer  behind  the  detec¬ 
tor  instead  of  a  finite  wall.  N3S-76  indicates  that  this  can  be  allowed 
for  by  deliberately  lowering  the  plotted  curves  for  small  X  ,  so  that 
the  resulting  exposure  rate  is  unity  for  X  -  0  and  H  =  3  feet  .  Thus, 
D2  ^  2W ( X , H )  (4) 

for  all  but  very  small  values  of  X  .  However,  there  are  some  minor 
approximations  inherent  in  the  above  formula  for  W  ,  so  that  Equation 
(4)  is  not  rigorously  correct,  even  for  very  large  values  of  X  .  (See 
NBS-42,  p.  3d' 

Situation  2-0  is  the  same  as  Situation  2,  except  that  the  mass  thick¬ 
ness  of  the  walls  is  zero.  The  detector  response  is  the  same  as  it  would 

be  in  the  open,  that  is, 
i 

°2-0  "  I  £,F.cose'  d(cos0)  =  L(K'  (5) 

J-l 

where  the  ^.-function  is  defined  in  NBS-42  ^Equation  27.2).  D 

<Ctmm  U 

obviously  equals  when  H  equals  three  feet. 

Situation  2-°°  is  also  similar  to  Situation  2,  except  that  the  mass 
thickness  of  the  walls  is  approaching  infinity.  This  is  to  be  regarded 
as  being  brought  about  by  maintaining  the  actual  thickness  but  letting 
the  density  approach  infinity.  in  comparison  with  the  value  of  D2,  D2_Q 
is  negligible.  However,  for  the  sake  of  comparison  with  exposure  rate 
for  other  configurations  having  walls  of  mass  thickness  approaching 
infinity,  one  must  consider  ir.  the  limit  the  variation,  with  respect  o 
direction,  of  the  differential  contribution  of  radiation  entering  the 
detector.  The  situation  is  one  which  is  extremely  difficult  to  establish 
experimentally,  but  there  are  theoretical  means  of  obtaining  radiation 
data  at  very  great  penetration  depths;  and,  since  the  angular  and  energy 
distribution  of  the  radiation  field  fends  to  seek  an  approximate  equilib¬ 
rium  at  great  depths,  there  is  reason  to  expect  that  a  description  of 
can  be  obtained  in  such  a  way.  The  primary  information  needed  is 

2-eo 
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the  directional  distribution  of  radiation  emitted  from  the  wall,  or  - 
looking  at  the  matter  from  the  detector  point  of  view  -  the  directional 
distribution  of  differential  exposure  contributions  at  the  detector. 

NBS-76  implicitly  assumes  that  the  directional  density  function  depicting 
the  radiation  emitted  from  the  walls  is  separable  into  components  depen¬ 
dent  on  vertical  (polar)  and  horizontal  (azimuthal)  angles,  respectively. 
Thus,  this  function  is  assumed  to  have  the  form: 

F(cos8,cp)  =  p (cp)  q(cos0)  ,  (6) 

where  0  and  cp  are  the  polar  and  azimuthal  angles,  respectively  rela¬ 
tive  to  a  frame  of  reference  with  a  vertical  polar  axis.  If  the  azimuthal 
angle  is  measured  relative  to  a  line  normal  to  a  wall,  then  NBS-76  gives 
as  an  expression  for  p(cp)  ,  based  upon  an  assumed  cosine  distribution, 


It  is  to  be  noted  that  the  adequacy  of  the  assumption  as  to  separability 
of  the  function  F(9,cp)  has  never  been  studied;  and,  even  if  it  is 
approximately  separable,  the  accuracy  of  the  cosine  assumption  for  p(cp) 
has  not  been  verified. 

Situation  3-0  involves  a  circular  structure  of  zero  mass  thickness; 
and  thus  D3_Q  is  the  exposure  rate  above  an  infinite  contaminated  plane, 
on  the  axis  of  a  decontaminated  circle  which  subtends  the  same  solid  angle 
at  the  detector  position  as  does  the  rectangular  base  of  the  structure 
under  consideration.  A  general  formula  for  the  exposure  rate  can  be 
expressed  as: 
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D 


3-0 


■t(H,cos0)  dfi 


(10) 


where  the  integration  is  over  the  solid  angle  subtended  at  the  detector 
by  the  walls.  Through  the  cylindrical  symmetry  of  the  situation,  this 
becomes: 


0 

D3-0  =  j 


COS0 

f 1 

<6(H,cos0)  d(cos0)  +  I  {.(H,cos0)  d(cos0)  (11) 


-cos© 


-cos© 


-t(H,cos@)  d(cos0) 


(12) 


This  formulation  is  not  exact  in  principle,  since  the  {.-function  is 
for  a  completely  contaminated  plane  and  a  collimated  detector,  while 
this  case  is  for  a  partially  cleared  plane  and  a  completely  uncollimated 
detector.  The  degree  of  approximation  is  quite  close,  however,  in  most 
practical  cases. 


Situation  3-°q  has  a  rectangular  plan  similar  to  that  of  the  structure 
itself,  but  it  is  considered  of  infinite  extent  in  the  vertical  direction. 
The  mass  thickness  approaches  infinity.  The  statements  relative  to  Situa¬ 
tion  2-oo  apply  in  general  here  also.  Under  the  same  assumptions,  the 
exposure  rate  is  found  to  be: 


j. 

L 


D  =  2  J  q(cos0)  d(cos0) 
3-oo 


/ 


TT 

2 


=  2(sin9  coscp 

8 


-i 


coscp  dcp  + 

0  *"0 
q(cos0)  d(cos0) 


<P« 


cos(p  dcp 


(13) 


(14) 


Situation  3!-oe  is  similar  to  Situation  3-»,  except  that  the  plan  is 


circular  instead  of  rectangular.  The  circle  has  the  same  area  as  that 
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specified  for  Situation  3-0.  NBS-76  makes  the  same  assumption  of  sepa¬ 
rability  of  vertical  and  horizontal  dependence,  and  in  this  case  an 
additional  assumption  is  made  that  the  vertical  dependence  is  similar 
to  that  of  the  skyshine  portion  of  the  ^-function,  that  is: 

q(cos0)  =  k  •  £(3* ,- ] cos0 1 )  ,  (15) 

where  k  is  an  appropriate  constant.  The  dependence  on  the  negative 
of  the  absolute  value  of  cos0  shews  that  the  contribution  from  above 
and  below  the  detector  plane  is  believed  to  be  symmetric  at  deep  pene¬ 
trations.  (One  should  note  that  here,  just  as  in  Situation  3-oo,  the 
wall  is  effectively  of  infinite  depth  as  well  as  infinite  height.)  For 
this  case  any  line  from  the  detector  is  perpendicular  to  the  wall,  and 
thus  from  Equation  (7),  p (cp)  equals  1/2  at  the  detector  position. 
Then 


2n 

D  s|  4. 

S'-00  I  2 

Jo 

0 

r 


2k  |  -t(3*  ,  cosS)  d(cos9) 


-1 


=  2nk  1  -M3',  cos0)  d(cos0) 

-1 

=  2nk  S(3 ' ) 


(16) 


(17) 

(18) 


where  S  is  a  function  defined  in  NBS-42  (Equation  27. 3^. 


Situation  4-0  differs  from  Situation  3-0  only  in  that  the  cleared 
area  within  the  "phantom''  structure  is  rectangular  rather  than  circular, 
oeing  similar  in  plan  to  the  structure  under  consideration.  A  general 
formula  for  the  exposure  rate  is  given  by 


D 


4-0 


£(H,cos0)  dQ 


(19) 


where  the  integration  is  over  the  solid  angle  subtended  at  the  detector 

(15  ) 

by  the  walls  of  the  structure.  The  technology  as  established  assumes 
the  identity  of  and  D4_q  ,  although  NBS-76  notes  that  in  principle 

this  is  not  true.  Calculation  of  D4_Q  can  be  made,  although  it  is 
tedious  because  the  limits  of  integration  are  not  simple  functions  of  the 
polar  coordinates  which  are  suitable  for  the  ^-function. 
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Situation  4-qq  is  the  situation  for  a  finite-sized  building  of  rectan¬ 
gular  plan.  It  is  the  same  shape  as  the  structure  under  investigation, 
but  it  has  walls  of  infinite  density.  The  angular  distribution  of  detec¬ 
tor  response  in  such  case  has  not  been  previously  defined,  but  in  principle 
it  can  be  obtained  by  theoretical  means,  if  one  can  show  that  deep  pene¬ 
trating  radiation  approaches  an  equilibrium  directional  distribution.  In 
terms  of  the  function  F(cos,cp)  introduced  in  the  discussion  of  Situa¬ 
tion  2-oo  t  the  exposure  rate  is  simply  given  by: 


D 


4-00 


F(cos9,cp)  dQ 


(20) 


where  the  integration  is  over  all  solid  angle  subtended  by  the  walls  at 
the  detector. 


Situation  4' -qq  is  similar  to  that  in  Situation  3'-oo  ,  except  that 
the  upper  and  lower  limits  on  the  wall  are  finite.  Thus  we  would  have, 
according  to  the  distribution  assumed  by  NBS-76: 


D  =  rrk 

4  "*oo 


{/ 

-C" 

0 

I  l(! 


t(3' ,  cos0)  d(cos0) 
COS0. 


3',cos0)  d(cos0) 


-cos© 


M- 


(21) 


=  TTk 


G(l-cos0.)  +  G(l-cos0  ) 
l  U 


(22) 


where  G  is  defined  in  NBS-76  (Section  4.4,  line  25).  The  validity  of 
this  equation  is  based  on  several  assumptions  as  noted  for  Situations 
3-oo  and  3'-»  . 

Situation  4  is  the  real  structure,  of  finite  dimensions,  actual  wall 
density,  and  rectangular  shape.  The  exposure  rate  is  expressible  as  an 
integral  over  all  directions  subtended  by  the  walls  at  the  detector  of 
the  directionally  differential  detector  response  function,  but  this 
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definition,  though  obvious,  is  not  directly  useful.  A  direct  and  exact 
calculation  is  not  easily  done;  and  it  would  require  an  elaborate  Monte 
Carlo  type  of  computation  not  suited  for  engineering  purposes. 


B.3  Discussion  and  Conclusions 

Since  R  is  defined  as  D^/D^  ,  the  value  of  R  may  be  established 
in  terms  of  the  exposure  rates  for  the  various  situations  already  dis¬ 
cussed,  according  to  the  following  analysis: 


One  may  consider  that  the  ratio  D  /D  ,  which  involves  true  wall  mass 

thicknesses,  can  be  interpolated  between  the  corresponding  ratios  for 

infinite  and  zero  mass  wall  thicknesses,  respectively,  that  is  , 

D,  /Dn  and  D  /D.  .  .  The  interpolation  factor  is  called  S  , 
4-00  2-oo  4-0  2-0  w 

and  thus 


D 


4 


D 


2 


NBS-76  assumes  that 


S  = 
w 


B  (X)-l 
P 

B  (X) 

P 


(24) 


(25) 


where  B  is  the  buildup  factor  for  cobalt-60  normally  incident  on  con- 
P 

Crete  (or  similar  material).  Such  assumption  involves  two  separate  hypo¬ 
theses:  (1)  the  only  significant  parameter  affecting  the  value  of 

is  the  actual  mass  thickness  of  the  structural  wails;  (2)  the  dependence 
is  given  by  the  ratio  of  scattered  to  uncollided  cobalt-60  radiation  which 
can  penetrate  the  walls. 


Thus,  one  obtains,  following  this  approach: 
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If  the  assumption  is  made  that 


D 

D 


4-00 

3~00 


4 '-00 


(28) 


which  involves  some  degree  of  approximation,  however  close,  then  Equation 
(1)  follows,  through  the  establishment  of  the  following  definitions: 


B 

w 


G 

s 


E 

s 


E 


d 


4'— CO 


00 


D 


3-0 


D 

D 

D 

D 

D 


2-0 

3~oo 

2- oo 

4-0 

3- 0 


(29) 

(30) 

(31) 

(32) 

(33) 


The  mathematical  descriptions  of  these  concepts  follow  from  the 

definitions.  Since  D,  is  unity,  B  is  obtained  by  the  same  formulation 

1  w 

as  D  (see  Equation  3  and  4  and  associated  text).  G  becomes,  from 

A  S 

Equations  (18)  and  (22), 


G 

s 


G(i-cos0.)  +  G(l-cos0  ) 

L _ _ _ _ u  —I 

2  S(3 ' ) 


(34) 


5 


l-cos0  ) 


G(l-cos0^) 


(35) 


since  S (3 ' )  is  approximately  equal  to  0. 1  .  From  Equations  (5)  and 
(12),  one  obtains: 
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G 


d 


£(H,cos0) 


L(H) 


d(cos0) 


(36) 


From  Equations  (9)  and  (14),  one  derives: 

E  =  sincp  +  cosep  .  (37) 

5  3  3 


From  Equations  (10)  and  (19),  we  see: 


E 


d 


£(H, cos0)  dn 


<L(H,cos0)  df) 


(38) 


where  the  integrations  are  over  the  solid  angle  subtended  at  the  detector 

by  the  walls  for  the  rectangular  and  equivalent  cylindrical  situations, 

as  indicated.  NBS-76  (Section  4.5)  suggests  that  this  ratio  can  be  obtained 

( 16  ) 

from  the  work  of  Hubbell  et  al.  .  Although  the  situations  considered 

by  Hubbell  differed  slightly  from  those  involved  here,  the  results  so 

obtained  would  probably  be  a  close  approximation.  Present  engineering 
(15) 

practice  ,  however,  is  to  Tssume  this  value  to  be  unity. 

It  is  to  be  noted  that  the  concepts  involving  essentially  infinitely 

thick  structures,  G  and  E  ,  involve  situations  which  are  probably 

s  s 

impossible  to  reproduce  experimentally.  Independent  determinations  of 
these  so-called  "non-wall  scattered"  factors^5\  or  even  their  product, 
by  experimental  means  seem  out  of  the  question.  Indirect  experimental 
determination  of  these  factors  by  measuring  the  overall  reduction  factor 
cannot  be  achieved  unless  all  the  other  factors  in  Equation  (1)  are 
accurately  known  and  the  various  assumptions  inherent  in  the  development 
of  that  formula  have  been  justified. 


The  important  assumptions  and  approximations  which  are  subject  to 
verification  or  possible  improvement  are: 

(a)  That  deeply  penetrating  radiation  can  be  split  into  factors 
depending  separately  upon  vertical  and  azimuthal  angles. 
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(b)  That  the  azimuthal  variation  for  deeply  penetrating  radiation 
follows  a  cosine  distribution. 

(c)  Tha;  the  vertical  variation  for  deeply  penetrating  radiation 
is  approximated  by  the  skyshine  portion  of  the  /-function. 

(d)  That  the  shape  factor  EJ  is  always  very  near  unity. 

d 

(e)  That  the  interpolation  factor  S  is  dependent  only  on  wall 

w 

mass  thickness  and  that  it  is  approximately  given  by  the 
formula  presented  as  Equation  (25). 

(f)  That  G  as  determined  on  the  basis  of  cylindrical  structure 

s 

situations  are  valid  for  rectangular  structures . 

One  might  ask  at  this  point  whether  it  is  worth  seriously  considering 
abandoning  the  whole  approach  to  ground  contribution  to  structure  pene¬ 
tration  typified  by  Equation  (1).  However,  it  must  be  observed  that  the 
present  formulation  has  been  shown  by  many  experiments  to  be  rather  suc¬ 
cessful  thus  far  and  is  now  deeply  ingrained  into  the  minds  of  many 

fallout  shelter  analysts.  The  edifice  thus  appears  to  require  bolstering 

* 

rather  than  replacement. 

B. 4  Recommendations 

It  is  recommended  that  il  further  work  on  ground  penetration  into 
blockhouses  be  undertaken,  it  be  directed  along  the  following  lines  in 
the  immediate  future: 

(a)  Make  further  efforts  to  pin  down  for  all  practical  thick¬ 

nesses,  either  experimentally  or  theoretically  or  both,  so  as  to  give 
confidence  in  definitive  results  drawn  from  such  work  within  an  uncer¬ 
tainty  on  the  order  of  perhaps  10?-  . 

(b)  Bo  further  work,  either  experimentally  or  theoretically  or  both, 

to  obtain  or  improve  values  of  G  and  E  for  phantom"  buildings 

a  (1 

(i.  e. ,  those  of  zero  wall  mass  thickness),  to  the  same  order  of  accuracy 
as  indicated  above. 

V.  Eisenhauer  disagrees  with  this  (personal  communication).  He  would 
prefer  a  completely  fresh  approach,  avoiding  the  interpolation  procedure. 
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(c)  Determine  E  and  G  (or  their  product  if  they  are  not  sepa- 

s  s 

rablv.’)  by  approved  theoretical  techniques,  to  within  the  same  order  of 
accuracy  as  above. 

(d)  After  the  above  factors  are  known  with  reasonable  accuracy, 
resume  experimental  blockhouse  work  to  look  for  the  best  values  of  S 

w 

and  to  determine  the  effect  of  various  parameter  changes  on  it.  This 

could  possibly  be  supplemented  with  Monte  Carlo  calculations  us'ng  the 

actual  blockhouse  configurations.  If  the  other  functions  are  known  within 

an  uncertainty  of  about  10%  ,  one  should  be  able  to  determine  J  to 

w 

within  about  20%  . 

It  is  not  implied  that  previous  experimental  work  in  these  matters 
has  been  wasted.  However,  it  is  believed  that  after  items  (a),  (b), 
and  (c)  have  been  accomplished,  previous  work  could  be  re- interpreted  in 
the  light  of  this  approach,  along  with  further  work  mentioned  in  sub- 
paragraph  (d). 

Some  justification  for  continuation  of  work  on  this  basic,  elemen¬ 
tary  problem  appears  to  be  necessary,  especially  in  view  of  the  apparant 
adequacy  of  the  approach  for  simple  geometries,  as  indicated  above.  Such 
justification  arises  most  strongly  from  the  fact  that  a  great  many  of 
the  problems  involving  more  complicated  type  of  structures  and  source 
distribut ion*  are  solved  by  methods  which  spring  from  this  basic  method 
and  assume  its  adequacy.  It  is  well  known  that  there  are  serious  dis¬ 
crepancies  between  experimental  values  and  engineering  predictions  in 
connection  with  many  of  these  complex  situations,  such  as  the  blockhouse 
basement  problem  t  in- and- down"',  the  protection  prediction  for  upper 
floors  of  multistory  structures  (’  i n-and-up*  ) .  the  mutual  shielding 
problem  (’’finite  source  plane"'*,  and  others.  In  order  to  solve  these 
difficult  problems  in  fallout  shielding  methodology,  tho  logical  frame¬ 
work  of  the  rjore  basic  part  of  the  technology  must  be  sufficiently  firm 
to  assure  that  the  troubles  in  more  complex  situations  do  not  stem,  even 
in  part,  from  these  shaky  foundations.  Therefore,  more  careful  attention 
in  research  should  be  devo’ed  to  strengthening  the  fundamentals  of  the 
present  fallout  shielding  technology. 
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